Astronomy & Astrophysics manuscript no. 16369 


©ESO 2013 


January 20, 2013 





The flare model for X-ray variability of NGC 4258 

T. Trzesniewski^, B. Czemy^, V. Karas^, T. Pechacek^, M. Dovciak^, R. Goosmann"^, and M. Nikolajuk^ 



7—i ' Institute of Physics, Jagiellonian University, Reymonta 4, P-30059 Krakow, Poland 
^ Copernicus Astronomical Center, Bartycka 18, P-00716 Warsaw, Poland 

' Astronomical Institute, Academy of Sciences, Bocnf II 1401, CZ-14131 Prague, Czech Republic 

Observatoire Astronomique de Strasbourg, F-67000 Strasbourg, France 
^ Faculty of Physics, University of Bialystok, Lipowa 41, P-15424 Bialystok, Poland 



in 

(N 



Oh 
I 

O 

> 
oo 



o 



X 



Received 20 December 2010; Accepted 18 April 2011 



ABSTRACT 



Aims. We study the variability mechanism of active galactic nuclei (AGN) within the framework of the flare model. We examine the case of 
Seyfert/LINER galaxy NGC 4258, which is observed at high inclination angle and exhibits rapid fluctuations in its X-ray light curve. 
Methods. We construct a model light curve based on the assumption of magnetic flares localized in the equatorial plane and orbiting with 
Keplerian speed at each given radius. We calculate the level of variability as a function of the inclination of an observer, taking into account all 
effects of general relativity near a rotating supermassive black hole. 

Results. The variability level is a monotonic function of the source inclination. It rises more rapidly for larger values of the black hole spin 
(Kerr parameter a) and for steeper emissivity (index j8 of the radial profile). We compare the expected level of variability for the viewing angle 
81.6 deg, as inferred for NGC 4258, with the case of moderate viewing angles of about 30 deg, which are typical of Seyfert type-1 galaxies. 
Conclusions. Highly inclined sources such as this one are particularly suitable to test the flare model because the orbital motion, Doppler 
boosting, and light bending are all expected to have maximum effect when the accretion disk is seen almost edge-on. The model is consistent 
with the NGC 4258 variability, where the obscuring material is thought to be localized mainly toward the equatorial plane rather than fonning 
a geometrically thick torus. Once the intrinsic timescales of the flare duration are determined with higher precision, this kind of highly inclined 
objects with a precisely known mass of the black hole can be used to set independent constraints on the spin parameter. 

Key words, accretion, accretion disks - galaxies: active - galaxies: Seyfert - X-rays: galaxies 



1. Introduction 

For many years the X-ray emission of active galactic nu- 
clei (AGN) has been known to be strongly variable (see e.g. 
iGaskell et al.l2006l: IUttlevl20O7l for reviews). However, the na- 
ture of this variability remains unknown. Most proposed mod- 
els have been based on a natural assumption that the emis- 
sion originates close to the accreting supermassive black hole 
and fluctuates on the dynamical timescale. Rapid var iability 
is am ong the arguments in favor of this interpretation (IKrolik 
199^ . The effects of general relativity are expect ed to play 



ies predict that non-stationary continuum models can dis- 
play relativistic effects also through the dependence of the 
level of variability on the inclin ation angle o f an observer 



(Zhang & Bao 1991; Abramowi cz etaLlI 19911: Karas 199 



an im portant role in shaping the observed signal (iKato et al 



1998) 



As the innermost part of accretion flow likely proceeds 
through some form of a disk, characterized by a (roughly) 
Keplerian profile of rotational velocity, the relativistic effects 
are expected to depend on the inclination angle of the ob- 
server with respect to the disk plane. These effects are only 
weakly seen in the stationary continuum models if the emis- 
sion is due to Comptonization, but they are very profound 
in the observe d shapes of spectr al features such as the iron 
Ka line (e.g. iFabian et alj 1 1995b . In addition, earlier stud- 



lFukuell2003l: ICzernv et al.H2004 . Therefore, the application of 
a specific model to both low and high inclination sources opens 
an additional possibility of model testing. 

In the present paper, we apply the flare model developed 
originally for the case of Seyfert galaxy MCG-6-30-15, which 
is a source at a moderate inclinatio n (Or, ~ 3 deg) with respect 
to the observer's line of sight (Fab ian et al ]|l 995, 2002), to the 
NGC 4258 galaxy viewed from the almost edge-on direction. 
This active galaxy is unique in several aspects: well-resolved 
maser emission from the nucleus allows an accurate mass de- 
termination of the black hole; the source is Compton thin (de- 
spite the high inclination), which allows us to measure the X- 
ray variability; rotation studies also show that the inner accre- 
tion disk follows Keplerian orbital motion very accurately. 

The paper is organized as follows. In Sect. |2] we summa- 
rize the general scenario of the flare model, introduce a con- 
venient parameterization, and adapt this scheme to the case of 
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NGC 4258. We discuss the main differences in the model set-up 
caused by our concentration in previous papers on low to mod- 
erate inclinations for the application to unobscured Seyfert 1 
AGN, whereas now we wish to apply the model to a highly- 
inclined source. In Sect. [3] we present our showing in particular 
the inclination dependence of the variability variance. We con- 
sider the role of the avalanche prescription, where the flares are 
mutually interconnected as they occur in families. Finally we 
summarize the results in Sect. |4]and present our conclusions. 

2. Flare model for X-ray variability of NGC 4258 

2.1. Description of the flare model 

The idea of X-ray emission coming from multiple loca- 
tions within/above accretion disks around a black hole is 
motivated by the important role of magnetic fields in the 
process of accre tion. T h e orig inal formulation of the flare 
model by Galee vTt al.l ( 1979 ) was followed by numerous 
papers (e.g. Pi Matteo 1998; Navakshin & Kazanas 2001' 



Merloni & Fabian ,2001; Collin et al.. ,2003; .Goosmann et al.l 
2007bh . A common theme to these models is an underlying 



assumption of a standard-type accretion flow driving the mag- 
netic field. A similar picture can also be developed in the ab- 
sence of a cold disk since in that case mul tiple shocks are ex- 
pected to form in t he hot accretion flow (e.g. lBoettcher & Liang 
1999HZvckill2003h . The attractiveness of the flare model is also 



supported by the close correlation between the radio and X-ray 
emission in radio quiet AGN, which is a phe nomenon typical 
of active stellar coronae dLaor & Behadl2008l) . 



In the present work , we generalize the model d eveloped 



by Czernv et al. ([2004) and Goosmann et al. ( 2006 ). We pa 



rameterize the distribution of the flares across the disk and 
the flare properties. We assume that the luminosity of a sin- 
gle flare following the initial onset decreases gradually with 
time. This differs from the previous paper where we adopted 
a rectangular profile for the flare light curve. The most im- 
portant modification is the introduction of a coupling be- 
tween the flar es through the presence of avalanches, a s dis- 



bv lPechacek et al.l ( l2008h . Thus, our model follows the general 



(IMineshige et al. 


1994 


1 Lvubarskii 1997; Kotov ef alJ 12001 


lArivalo & Uttlevll2006 


; Titarchuk & Shaposhnikovl2008b. 



Spontaneous flares originate at a certain outer radius, ^spont- 
Its value can be rather distant from the center, as this seems 
to follow from the discussion of well-studied cases o f major 
flares in AGN (e.g. lPonti et al.ll2004 l; iGoosmann et"aDl2007a). 
On a more theoretical side. TUzdensk v & Goodm an (2008) de- 
scribed a specific scenario of a magnetized corona above a 
turbulent accretion disk. Magnetic loops rise above the disk, 
where they become sheared by differential rotation of the field 
line foot-points until reconnection occurs. These events are re- 
sponsible for the occurrence of primary flares and the subse- 
quent dissipation of magnetic energy. 

Spontaneous flares are accompanied by avalanches of sec- 
ondary flares developing subsequently inside Rspom- The as- 
sumptions about the parent flare localization and the grad- 



ual avalanche progression towards smaller radii go beyond the 
standard avalanche model, which does not invoke any particu- 
lar form of the spatial distribution. General arguments, based 
mainly on the diffusion equation for the propagation of the 
events, suggest that the properties of secondary flares are de- 
scribed by power-law distributions. In this way, the onset of 
variability occurs at large radii and induces secondary fluctu- 
ations further down the accretion flow, in the inner regions of 
the disk. 

In our paper, the flares are assumed to vanish at the in- 
ner radius of the computational zone, Rin, which corresponds 
to the final inflow of perturbations below the innermost sta- 
ble circular orbit (ISCO), and eventually below the black hole 
horizon. The assumption about the exponential decay of indi- 
vid ual flares is another improvement of our previous scheme 
(see ICzernv et al. I l2004 where we assumed a common charac- 
teristic lifetime for the duration of all flares which then switch 
off abruptly). Onset of flares followed by their rapid but gradual 
decay may be analogous to processes in solar flares, where sim- 
ilar mechanisms of heating vi a magnetic reconnection operate 
(e.g.. lAschwanden et al.lll998l) . 



2.2. Parameterization of the model 

We briefly summarize the adopted parameterization. Flare 
events are characterized by radius from the black hole and the 
corresponding time of the occurrence. During their lifetime, the 
flares follow a coplanar Keplerian motion of the accretion disk. 

As we wish to take general relativity effects into account, 
the flare event coordinates need to be defined in an appropri- 
ate manner. For the purposes of our investigation, the Boyer- 
Lindquist (f, r) coordinates suit well, because all flares are as- 
sumed to be distributed above the ISCO, so these coordinates 
are non-singular over the whole computational domain. An in- 
dividual flare is located at the radius r and characterized by the 
moment of its birth, fbkth. The flares rise instantaneously and 
decay with time t exponentially (with a sharp final cut-off) 



cu ssed e.g. by Poutanen & Fabian ( 19991) . Zvcki (2002), and - 



fi(r) exp (-e ^-p^) , fbiith < t < tunh + T(r), 



(1) 



0, 



t > fbiith + T(r) or f < 0. 



The amplitude /i(r) of a flare depends on the flare location, and 
its dependence on the radial coordinate has the form 



\'-o 



-P 



(2) 



where /o is a normalization constant and /? is a model parame- 
ter. The lifetime of the flare is assumed to be related to the flare 
location as 



(3) 



where tq and 6 are the model parameters. The scaling radius ro 
was set at 187?g. 

In our picture, we assume that each flare can lead to a new 
flare with the probability given by a Poissonian distribution 
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around a mean value yu. Each new flare is always generated 
closer to the black hole than the parent flare: the location of 
the new flare, rnew, is derived from the radii between the parent 
flare radius, r, and the radius ir (0 < 5 < 1), where the radial 
distribution of the probability is defined as 



P('-newk) = n{)(r)r^ew, 



(4) 



where no(r) is the local normalization constant depending on 
the position of the parent flare, and y is a model parameter The 
value 7 = describes the uniform distribution across the disk 
surface. The time of birth of a new flare is drawn uniformly 
from the time bin between the birth of the parent flare, t, and 
/deiayf, where the dimensionless factor /delay is of the order of a 
few units. 

The level of variability is set by global parameters: the 
number of spontaneous flares, A^f, and the total duration of 
the light curve, Ttot. The size of the single emitting region 
is then identical for all flares. It is fixed by the number of 
flares present in the source and the average luminosity L of the 
source. Flares are in Keplerian motion around the black hole. 
All general relativity effects connected with light propagation 
toward an observer are calcu lated using the ky code (for de- 
tails see Dovciak et al. ]L2004a"b'). For that purpose, all emitting 
regions are assumed to be projected onto the equatorial plane. 

The shape of the locally emitted spectrum is parameterized 
as a power law of a given slope, and the reflection component 
with arbitrary normaUzation is usually inc luded in the model 
dCzemv et alj|2004 iGoosmann et al.ll2006h . 



2 galaxy because the emission lines seen in the polarized light 
have widths appropriate for this class of objects. 

The source was classified in X-rays as an obscured 
low-l uminosity active nucleus by the ASCA observa- 



tion dMakishimal Il994 . It was observed at X-ray wave- 
lengths by two high-spatial resolution instruments: Chandra 
dYoung & Wilson 2004) and XMM-Newton (Pietsch & Read, 
2OO2I) . More recently, results derived from Suzaku ob- 



2009 



serva tions ha v e also been reported (iRevnolds et al. 
Yamada et al. 1 l2009|). The nuclear emission is s t rongly 



variable (jMarkowitz & UttlevI 120051: IRevnolds et al.l 12009 



Yamada et al.ii2009i) . It follows that the surrounding medium 



cannot be Compton-thick along the line of sight from the en- 
tire X-ray producing volume. The obscuration instead causes 
a partial ab sorption with the absorbing column ~ lO^-' cm~^ 
dYang et a l. 2007; Yamada et al. 2009). The obscuring material 
can be localized mainly towards the equatorial plane rather than 
form a geometrically thick torus. 

The slope of the X-ray spectrum is typical of Seyfert galax- 
ies, where the value F ~ 1.7-2.0 depends on the model details 



(iPietsch & ReadI 120021: lYang et al.1 120070 . The cold and nar- 
row iron line comes from the outer region, whereas no broad 
line wa s required to fit the XMM-Newton spectrum of this 



Yamada et al 



source (lYang et al.ll2007t IRevnolds et all 120091 
(2009j) confirm a significant X-ray variabiUty and conclude 
that a classical geometrically thick torus does not seem to be 
present in this object, but that instead, the obscuration may take 
place along the line of sig ht clos e t o the e quatorial plane (e.g., 
iBao & Stuchlfklll992l:iKaras & Baoill992l) . 



2.3. Properties of NGC 4258 

NGC 4258 (M 106) is a barred spiral galaxy with a low- 
luminosity, type 1.9 Seyfert nucleus, which is also classified 
as a LINErH The nuclear source has a luminosity of L ^ 
lO-'^LEdd (iFruscione et al.ll2005l) and L2-10 = 3.31 x 10"*° erg/s 
in the 2-10 keV band. The geometrical distance to NGC 4258 
can be inferred from the orbital motion of H oO masers in its 
nucleus, as demonstrated for the first time by Herrnstein et al 
19991) . obtaining (7.2 + 0.3) Mpc. The mass measurement is 



reliable because the maser rotati on curve is most accurat ely de- 
scribed by the Keplerian profile (iHermstein et al 1 l2005h . since 
the radius of the inner masers, measured with the VLBI, cor- 
responds to 2.8 mas, high-velocity masers yield a mass f or the 
centi-al black hole of M = 3.78 x IO^Mq (lMartinll2008l) . The 
corresponding Eddington luminosity is Lead = 47.5 x 10"*"* — , 
hence L * 0.00475 x 10^"* erg/s. 

The masers indicate that the angle of inclination of the inner 
disk is high with respect to the observer line of sight: 0„ ~ 
81.6°, where 90° rep resents edge-on orientation ( Martinl2008l) . 
IWilkes et al.l ( Il995h find that optical lines are strongly linearly 
polarized and that the position angle coincides with the plane of 
the maser disk. Intrinsically, the source is likely to be a Seyfert 



1 See The NASA/IPAC Extragalactic Database (NED), which is 
operated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and Space 
Administration. 



2.4. X-ray excess variance in NGC 4258 and tine 
enhancement ratio 

The X-ray flux of NGC 4258 is highly variable. The level of 
this variability is con veniently characteriz ed by the dimension- 
less excess variance (Nandraetal ] ll997ah 



1 ^ 

(rms/mean)2 = — J^^Xi - - uj]. 



(5) 



where is the number of (equally spaced in time) data points, 
Xi is the count rate in a given time bin, jj is an unweighted 
arithmetic mean of Xi, and cr, is the measurement error (pre- 
dominantly consisting of Poissonian noise). The error level de- 
pends on the duration of the light curve. This dependence is a 
typical property of the so urces dominated by th e red noise, as 
discussed for example bv lVaughan et al.l (l2003h . 

The source was observed sev eral times with the A SCA 
satellite. These data were used bv iNikolajuk et al. I (l2009h . who 
obtained 



(rins/inean)l_, = 6.215!^:|^xl0 



(6) 



for several light curves of typical total duration T - 31 300 sec- 
onds, consisting of 50 bins and the duration time of a single bin 
fb - 626 s. 

Since the black hole mass in NGC 4258 is known, this value 
can be compared to the expected variance from the relation be- 
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tween the black hole mass and the X-ray excess variance, that i s 
characteristic of Seyfert galaxies (INikolaiuk et al ]l2006ll2009l) 



Furthermore, 



(rms/mean) - C 



M ' 



(7) 



where C= 1.92M0S"'. 

For the black hole mass 3.78 x 10^ Mq (see Sect.|231l and 
the light curve parameters T and fb as above, we obtain 



(rms/meantp^^,,, = 1.52x10-1 



We thus define the enhancement ratio 



•7? = 



(rms/mean)\ , 

^ ' ^observed 



(rms/mean)^ , , ' 

^ ' ^expected 

and in the case of our source this ratio is equal to 
7? = 4+l. 



(8) 



(9) 



(10) 



This (rather large) error arises mainly from the systematic error 
in the constant C in Eq. 12.41 

We expect that the enhancement is due to the inclination 
of the source, which in the case of NGC 4258 is significantly 
larger than the average inclination (about 30 deg) of sources 
usual ly used to derive the Eq. (12.4b (see e.g. .Nikolaiuk et al. 
2006h . 



2.5. Adaptation of the model to NGC 4258 

In our calculations, we used values of the black hole mass and 
an AGN luminosity appropriate for NGC 4258 given in the pre- 
vious subsection. For the physical structure of the model, be- 
cause of the lack of hard evidence of a disk in the vicinity of 
a black hole in the NGC 4258 nucleus, we assumed that only 
primary emission, from the flares, is present (no cold accretion 
disk, therefore, no spots are induced by flares on the disk sur- 
face). We assume that this emission is concentrated all around 
the equatorial plane. 

Some model parameters were constrained to obtain a pro- 
file of the X-ray power-spectral density that is close to that 
for MCG-6-15-30 as a general representation of the observed 
shape for AGN. The values of the model parameters that we 
did not vary are as follows: 5 = 1.5 (assuming that flare dura- 
tion scales with local Keplerian timescale), y = 1.5, = 0.2, 
/delay = 7, // = 1.1, and 7?out - 200. Dimensionless dis- 
tances are given in units of the black hole gravitational radius, 
Rg = GM/c'. Our model of course allows us to change the 
above-mentioned parameters if the physical reasoning requires 
us to use different values. We tested the model dependence on 
the luminosity concentration given by jS (see Eq.|2]l, on the pro- 
portionality constant of the flare duration, tq (see Eq. O, and 
on the Kerr parameter, a. 

The quantities in physical units, namely, the mass, spin, and 
distances, are related to the corresponding values in geometri- 
cal units 



M 



^phys 1.477 X lO-'^cm' 



flP"^^ = ca. 



(11) 



a 
M 



phys 



r 

M 



rPhys 



(12) 



We employ geometrical units in numerical simulations, but 
convert them to physical units interpreting the results. To obtain 
frequency in physical units [Hz], the relation is af^^'^ - cu). The 
geometrized frequencies scale as M"'. Their numerical values 
must thus be multiplied by a factor 



c A ( M , 
=(3.231x10^) [Hz]. 



(13) 



To calculate the light curves, we followed the ASCA setup, as 
the source was observed several times with this satellite. This 
helps us t o determine the excess variance to a sufficiently high 
precision ( Niko laiuk et al .l2009h . Therefore, in our calculations 
we applied the total duration of the observation T - 31 300 sec- 
onds with 50 bins and the duration time of a single bin th = 
626 s. 

Since the X-ray spectrum of NGC 4258 does not show the 
reflection component, the reflection must be either weak or ab- 
sent. Therefore, in the present study we simplify the model and 
neglect the reflection. In this case, the additional parameter re- 
lated to the normalization of the reflected component is unnec- 
essary. The results in this case even more remarkably do not 
depend on the assumed slope of the power law emission. 

For all model light curves, we calculated the variance and 
then tested its dependence on the inclination angle, while vary- 
ing the Kerr parameter of the black hole and the radial distribu- 
tion of the flare luminosity, The latter determines whether the 
emissivity is more or less concentrated toward the black hole. 

3. Results: Dependence of variance on inclination 

We use the aforementioned model to calculate the normalized 
X-ray excess variance, (rms/mean)^, paying special attention 
to the relation between this quantity and the observer inclina- 
tion. The flare model variability is very sensitive to the incli- 
nation, so we first explore this dependence while several other 
model parameters are kept fixed. In particular, we set the black 
hole mass to the appropriate value for NGC 4258. As for the 
black hole rotation, since we have no constraints for the spin 
parameter in NGC 4258, we also keep it free and analyze the 
level of variability for different values of a. 

In Figure [T] we plot the X-ray variance for three values 
of y6 assuming a rapidly rotating black hole (Kerr parameter 
a = 0.95). The value of jS ^ 1.5 corresponds to a shallow 
emissivity profile. The variabihty level increases with the in- 
clination, first slowly and later sharply. The rise is the result 
of the beamed Doppler boosting at high inclinations, which 



is well-known from previous papers (e.g. ICzernv et al.l 12004 



Dovciak et alj |2008). The rise is faster for larger values of j8, 
which implies that there is a larger contribution to the total lu- 
minosity of the source arising in the inner region. As a conse- 
quence, the relativistic effects are strongly enhanced. 

The aforementioned increase in the variance with inchna- 
tion is slower and the overall variability level is lower if the 
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Fig. 1. The dependence of the normalized variance on the in- 
clination angle of an observer for three emissivity profiles (3 for 
the Kerr black hole a - 0.95. The vertical line indicates the 
inclination of do - 81.6°. 




; [degj 



Fig. 2. The dependence of the normalized variance on the incli- 
nation angle for the same values of /? as in the previous figure, 
but with a non-rotating (Schwarzschild, a-Q) black hole. 



black hole is not rotating. In Figure |2] we show the inclina- 
tion dependence of the normalized variance for the same set of 
radial profiles (3 as in Fig. [1] but now with the Schwarzschild 
black hole. The ratio of the variances for j8 = 3.5 at 0o - 81.6° 
and 0o - 30° is equal to 3.5 for an a = black hole, but in- 
creases to 3.7 for an fl = 0.95 Kerr black hole in the presented 
plots. 

The rise in the variance with the inclination angle is not 
a result of the statistical errors in simulations. It is well 
known that a single simulation of short timescale is strongly 
affected by the pow er leaking from the lower frequencies 
jVaughan et al ]|2003h . lowever, in Figs.[T]and|2]the sequences 
for fixed values of fS were calculated for the same realizations of 
the same statistical distribution. The rise in variance is entirely 
due to the change in the viewing angle. If we use different ran- 
dom realizations of the flare distribution the trend is less clear 
because of large statistical dispersion in the adopted length of 



the light curve. A single variance (in simulations as well as in 
the actual data) is determined with an accuracy of a factor of 
two. 

Therefore, to show the model predictions to higher accu- 
racy, we had to extend the simulated light curves by a factor of 
eight. For these longer light curves, we calculated the variance 
enhancement, defined as the ratio of the variance seen at 8 1 .6° 
to that at 30°, as a function of the Kerr parameter The result is 
shown in Fig. [3] for three values of the flare time-scale duration. 
We also plot the ratio of the observed variance, given in Eq. (|6|l, 
to the variance expected from Eq. (4) of lNikolaiuk et al.l(l2009l) . 
This ratio is found to equal 4.0, its standard deviation (Icr) er- 
ror coming from the errors in both the variance and in scaling 
constant. Because of a lack of information about the black hole 
spin, the observational constraint is indicated by straight lines 
irrespective of the value of a. 

The variance enhancement, in general, increases with the 
Kerr parameter but the detailed profile depends on the assumed 
scaling factor of the flare lifetime. The variability is yet greater 
when the flare lifetime close to the black hole is comparable to 
the time spent by the flare in the region of the highest Doppler 
boosting. The angular extension of this region ( as a fraction of 



2n) decreases as the Kerr parameter increases (IPovciak et al. 
I2004al) . 

The enhanced variability is consistent with the expectations 
of the relativistic enhancement, within the framework of the 
flare model. A flare duration of longer timescale, tq = 10^ s, is 
indicative of a non-rotating black hole, while the flares scaled 
down to T() = 10^ s probably correspond to a spinning black 
hole with a > 0.5. However, taking into account the large error 
in the observed variance and the lack of a priori knowledge 
of the flare timescale we cannot at this stage make any firm 
conclusion about the black hole rotation in NGC 4258. 



4. Discussion 

In the flare model of AGN X-ray variability the change in the 
X-ray flux is caused by a combination of two main effects. The 
first is the intrinsic variability of both a single flare and the flare 
distribution. The second is the variation caused by the relativis- 
tic effects of the flare orbital (Keplerian) motion. Analyzing the 
change in the level of variability with the inclination we can 
disentangle those two effects. 

We chose NGC 4258 as a target of our study because of 
its exceptional properties. The source is highly inclined, seen 
almost edge-on, but is nevertheless Compton thin and the vari- 
able primary emission is still transmitted effectively through 
the material located at the source equatorial plane, although the 
X-ray emission is significantly absorbed. In addition, the mass 
of the central black hole in this object is measured accurately 
thanks to the spatially resolved water maser emission. 

On the other ha nd, we modified our original flare model 
(ICzernv et al.ll2004l) by including the exponential decay of in- 
dividual flares and the effect of an avalanche-type development 
of flares. On the other hand, we neglected the cold disk reflec- 
tion. In NGC 4258, any manifestations of reflection features, 
such as the broad relativistic iron line and temporary spots on 
the disk surface, are strongly suppressed. We thus understand 
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Fig. 3. The ratio H of normalized variance at the observer in- 
clination angle 0o - 81.6° with respect to the case 6o - 30°. 
The dependence on the dimensionless Kerr black hole param- 
eter a is plotted for the assumed emissivity profile fi - 3.5, 
for three values of the flare lifetime scale fo, as indicated in the 
plot. Three horizontal lines give the enhancement ratio H for 
NGC 4258 given in Eq. ( fTOl l (middle line), and the range of er- 
rors (upper and lower lines). The simulated light curves were 
longer by factor 8 compared to those used in Figs. [T}j2] 

the situation by presuming that the central corona is the source 
of enhanced variability and the place where flares originate. 

After developing a refined version of the model, we ana- 
lyzed the dependence of the X-ray variance on the inclination 
angle. We noticed that the rise in the variability with the in- 
clination angle is faster if the emissivity is more concentrated 
towards the black hole, which is the case for larger values of 
the j6 parameter We compared the expected level of variability 
for the viewing angle 81.6 deg, as inferred for NGC 4258, with 
the case of moderate viewin g angles 6„ ^ 30 deg , which are 



typical of Seyfert 1 galaxies (iNandra et al.lll997bi |2007|) . Our 
model predicts that the variability increase is monotonic with 
inclination. 

The observed variability level in NGC 4258 is indeed 
higher by a factor consistent with the measured value of the 
central black hole mass, which was determined with high ac- 
curacy in this source from water maser emission. The measured 
X-ray excess in this so urce is equal to cr- - 6.215^2 350 ^ 10 
jNikolaiuk e t al.ll2009l). The expected variance from the Eq. (4) 
of iNikotaiuk et alj ( 20091) is a factor of four lower (the afore- 
mentioned formula is satisfactory for typical low-inclination 
Seyfert galaxies). This means that the flare model of variabil- 
ity would be consistent with this object provided that the X-ray 
variance for a highly inclined source is a factor of four higher 
than it is for low-inclination sources for which the relation was 
scaled. 

Several parameters determine the outcoming light curve in 
our model. The interplay between the parameters is in general 
quite complex because the signal depends on the initial distri- 
bution of parent flares, their subsequent propagation across the 
accretion disk towards the inner edge, as well as the avalanche 
mechanism generating the secondary flares. To capture this en- 



tire process, we introduced the flare rise time fbhth in eq. ([TJ, 
the index of radial distribution of flare amplitudes fi in eq. (|2]l, 
the radial profile of flare lifetime S in eq. Q, and the index of 
radial probability distribution y in eq. (HJi. Furthermore, the sec- 
ondary flare generation is described by the parameters s, /delay, 
fi, Rin, and /?out, which define the development and the gradual 
decay of the avalanches. 

Despite the complicated picture given by the number of free 
parameters described above, it is interesting to note that only 
in special (but quite representative) cases do some parameters 
play an important role. This is in particular the case of a sta- 
tionary distribution that is expected to arise when the secondary 
flares occur over the whole range of radii in the accretion disk 
(5 —> 0). One can then check that the stationary situation de- 
pends sensitively on fi (i.e., the mean value of secondary flares) 
but the degree of the dependence on other parameters is much 
weaker. 

5. Conclusions 

The above-mentioned rise of variability agrees with the model 
prediction for a steep emissivity profile and the flare nor- 
malization timescale of 10^-10"* s. Higher values of the flare 
timescales are consistent with the low values of the black hole 
spin in NGC 4258, whereas shorter timescales are consistent 
with a rapidly spinning black hole. 

The observed enhancement in variability supports the view 
that the X-ray emission is generated close to the black hole and 
subject to strong general relativity eff'ects. This does not im- 
ply that the emission originates preferentially from magnetic 
flares above the disk, or that it instead comes from individual 
shocks in an optically thin inner flow that undergoes Keplerian 
motion. In the absence of a reflection component, our model is 
unable to diff'erentiate between the two aforementioned scenar- 
ios. However, if the inflow is almost spheri cal, as in the mod- 
els of Shrader & Titarchukl(ll998.) and .lshibashi & Courvoisier 



(l2009h . the variability enhancement is not expected for high 
inclination sources. The observed enhancement of variability 
of NGC 4258 thus implies that there is a significant azimuthal 
motion in the X-ray emitting plasma of the inner accretion disk. 

Once the intrinsic timescales of the flare duration are deter- 
mined with higher precision, the highly inclined objects with 
precisely known black hole masses can be used to set inde- 
pendent constraints on the spin parameter The highly inclined 
AGNs that are not Compton thick and at the same time ex- 
hibit maser sources, such as NGC 4258, do exist (even though 
they are rather rare) and should be ideal candidates to apply the 
method described in this paper We will however need to obtain 
a better coverage of their varia bility properties. It app ears that 
unabsorbed Seyfert 2 galaxies (Brightman & Nandra '20081 or 
even dust-poor quasars ( Hao et al.. 2010) form a category of 
sources that could be suitable candidates. 
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